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Abstract. We propose a new cluster deprojection algorithm for recover- 
ing cluster structure along the line of sight (los) based on the Richardson- 
Lucy (RL) algorithm, which recovers a nonnegative theoretical distribu- 
tion from a given projection. To optimize our reconstruction we combine 
X-ray and Sunyaev-Zeldovich (SZ) data with gravitational lensing maps. 
The combination of these three data types constrains the structure of 
rich clusters of galaxies. To validate our new approach, we apply the 
algorithm to gas-dynamical cluster simulations. A reconstruction using 
the lensing potential is shown. 

1. Richardson-Lucy deprojection algorithm 

Resolving the los-structure of 3d objects can be formulated as inversion problem, 



where (j)(C) is the function of interest, ip(x) is the function accessible through 
measurement, and the integral kernel P(x\C) is normalized and non-negative as 
i/j(x) and (p(C) represent probability distribution functions. 

Using Bayes theorem for conditional probabilities, a two step iterative re- 
construction algorithm (Lucy, 1974; Lucy, 1994) 



for approximating 0(£) can be derived, with <jf +1 being the (r + l)-st estimate. 
For systems with axial symmetry inclined with respect to the los an explicit 
expression for P{x\C < ) was derived by Binney et al. (1990). 

2. The proposed algorithm 

The three observables depend differently on the los-structure of the 3d gravita- 
tional potential (Reblinsky & Bartelmann (1999)): 
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Figure 1. Reconstruction of the 3d gravitational potential (dis- 
played is —(f)) from the lensing potential. Inclination angle i = 30°. 
Left: the result after 7 iteration cycles; Right: the original cluster 
potential used to compute the input data, i.e. the lensing potential tp. 
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For the X-ray case and the SZ-effect an isothermal gas distribution is as- 
sumed. The inclination angle i and the temperature T are input parameters. 
This leads to the following combined two step reconstruction formula: 



Step 1: F ni = J f ni (x,y)P(x,y\R,Z)dxdy 

Step 2: 4>r+i oc 
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where the f ni are defined as f ni = ^-,f n2 = ^-, and f ns = . 



3. Application to test data from gas dynamical cluster simulations 

Using the lensing potential ij)(x,y) as input data, the characteristic features of 
the true 3d potential (f>(R, Z) are recovered very well after 7-8 iteration cycles. 
We also find that the reconstruction is insensitive to the initial guess, even 
(j)(R, Z) = const can be used as in Fig. 1. 
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